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ABSTRACT 

' Context. Direct detection of a planet around a star and its characterization for identification of bio- tracers in the mid-IR, requires 

a nulling interferometer. Such an instrument must be efficient in a large wavelength domain in order to have the capability to detect 
O ! simultaneously the infrared spectral features of several bio-tracers: C0 2 , 3 and H 2 0. 

Aims. A broad wavelength range can be effective, provided that an achromatic phase shift of n can be implemented, with an accuracy 
good enough for achieving a deep nulling at all considered wavelengths. A new concept for designing such an achromatic phase 
+3 \ shifter is presented here. The major interest of this solution is that it allows a simple design, with essentially one device per beam. 

Methods. The heart of the system consists in two cellular mirrors where each cell has a thickness which introduces, for a given central 
wavelength, a phase shift of (2k + 1)tt or of 2kn on the fraction of the wave it reflects. Each mirror is put in one of the collimated beams 
of the interferometer. Because of the odd/even distribution, a destructive interference is produced on axis for the central wavelength 
I t when recombining the two beams. Now, if the number of cells of a given thickness follows a rather simple law, based on the Pascal's 

triangle, then we show that the nulling is also efficient for a wavelength not too far from the central wavelength. 
^-l. ■ Results. The effect of achromatization is the more efficient the larger the number of cells is. For instance, with two mirrors of 64x64 

cells, where the cells phase shift ranges between -6n and +6n, one reaches a nulling of 10~ 6 on a wavelength range [0.6/t , 1.25/t ], 
i.e. on more than one complete octave. This is why we claim that this device produces a quasi- achromatic phase shift : especially, it 
could satisfy the specifications of space mission as DARWIN. In a second step, we study the optimum way to distribute the cells in 
the plane of the pupil. The most important criterion is the isolation of the planet image from the residual image of the star. Several 
algorithms are presented, one being especially efficient and we present the nulling performances of those various configurations. 

O 1 

Key words. Instrumentation: interferometers - Techniques: high angular resolution - Techniques: interferometric - Space vehicles: 
\ instruments - (Stars: ) planetary systems 
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1. Introduction unity, then it can be detected - in principle - in a much more 

$-H ' ,_ 10 efficient way, especially because of the reduced photon noise. 

' There is a huge contrast (10° iU ) between a star and a planet that 
orbits it, so that directly detecting the planet requires to cancel as Obtaining an achromatic n phase shift is one of the main keys 

much as possible the stellar light. Today, it is generally admitted of success because a) the wavelength domain is broad where 

that there are two main instrumental paths to achi eve this goal the set of a reputed unambiguous spectroscopic signature of life, 

: one is coronagraphy on a single telescope (e.g. lGuvonll2007l) namely C0 2 (15 /mi), O 3 (9.5 um) and H 2 (5-8 //m), is to be 

- generally in the visible to near-IR domain -, the other one found (Mger et al.||l996|; |Ollivier||2007|), b) in photons starving 

is nulling interferometry - proposed essentially in the thermal experiments such those considered here, a design based on a 

infrared - using at least two telescopes recombine d in a clever monochromatic detection is out of question, 
way (e.g. | Leger et aLl Eg^ |Woolf & Angelll997|). This paper Several so i utions have been presented in order to approach 

pertains to the second category. | Bracewell| 4l978|) was the first to an ac hromatic n phase sh ift in a large domain of wavelength : 

propose the concept of a nulling interferometer, where the light for a review? see |Rabbial d200l . The three main solutions are 

collected by two telescopes is coherently recombined, after that the s tacking of dielectric p lates of different indexe s and disper- 



a Ti phase shift be applied on one of the two arms of the interfer- s ion dMieremet et al J2 000). the crossing of a focus dRabbia et al.1 
ometer. The system of fringes projected onto the sky shows then [2005j) and th e use of a periscope opti cal train to reverse the 
a central dark fringe : if the star image is put on this central dark electric field dSerabvn & Colavital [2001). The first solution pro- 
fringe, it is strongly attenuated (actually, the stellar photons are vides only a qua si-achromatization since it relies on an approx- 
all sent m the second output of the interferometer). If, at the same imate linearization of the optical path vs wavelength. The two 
time, the planet is on a bright fringe, i.e. a region of transmission last metnods provide in principle an intrinsic achromatic phase 

shift, but which is strictly limited to n\ instrumentral effect (dif- 

Send offprint requests to: D. Rouan e-mail: daniel . rouan@obspm.fr ferential polarization for instance) may also affect the resulting 
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null. The Two first methods have the characteristic of introduc- 
ing a specific sub-system in one arm only of the interferometer, 
resulting in an assymetric design with the drawback of a diffi- 
cult balancing of the two arms. Finally one notes that all those 
solutions require the use of several optical components, requir- 
ing in general to fold the beams several times. Is there a way of 
reaching a n phase shift with a simple optical component, while 
keeping the interferometer design symmetric. We present in the 
following a n ew conc e pt tha t has th is ambi ti on. 

Recently iRouanl (120031 12004 I2006L l2007h presented a 
method to define interferometric configurations of telescopes 
that provide a n nulling function for any even value of n. 
Achieving a rather large value of n is important in order to pro- 
duce a trough shape of the nulling function and thus a good can- 
cellation of the leaks due to the finite angular size of the star. 
The principle i s based on a remarkable property discovered by 
iProuhetl (Il85ll) about a peculiar partition in two sets of the N 
= 2 L first integers, done according to the Prouhet-Thue-Morse 
sequence. This property is such that, selecting half of the N tele- 
scopes according to this sequence and applying to the output of 
each one a n phase shift, makes possible the cancellation of the 
(L-l) first terms of the developement in 6 of the recombined am- 
plitude, thus leading to the wished 6 2L behaviour of the intensity. 

Here, we remain in the same spirit while tackling the ques- 
tion of the achromatisation of the n phase shift. The idea is to try 
to apply a similar principle of cancelation of the first terms of 
a Taylor's development of the amplitude, but now with respect 
to the relative wavelength difference AA/Aq. The central com- 
ponent of the new concept is a chessboard mirror, i.e. a mirror 
made with a pattern of cells of different thickness. 

In this first paper, we consider only the case of an inter- 
ferometer with two square pupils telescopes and we establish 
mainly the theoretical grounds on which this new concept lies. 
We show that there is indeed a peculiar distribution of the cells 
thickness on the chessboard mirror that fulfils our goals of a 
quasi-achromatization of the n phase shift. We also show that 
an adequate distribution of the cells on the mirror's surface can 
add more nulling power. We depict several configurations, espe- 
cially one, that have good properties with respect to the detection 
of a planet. We then study the performances that a perfect device 
would feature and show that it would satisfy the specification 
of the currently propos ed missions, such as Darwin for instance 
(iLeger & Herbstll2007h . Several appendixes detail the formalism 
and the mathematical background for the readers who will be 
keen on a deeper investigation. A second paper will address the 
questions of the performances in a real world where a device 
cannot be perfect, and will present the first results of an experi- 
mental demonstrator. 



2. The concept of the chessboard phase shifter 

2. 1 . Formulation of the problem 

Let's consider an interferometer made with two identical afocal 
telescopes whose output beams are recombined in the so-called 
co-axial mode thanks to a beam-splitter. The cartoon of FigQ]de- 
picts such an optical scheme. To simplify the design and forth- 
coming computations, we assume that the shape of each tele- 
scope's pupils is limited by a square. In a collimated beam is- 
sued from the first telescope, we introduce a square cellular mir- 
rors (alternatively transparent plate could be used). This cellular 
mirror has a chessboard design in the sense that it features nxn 
square cells, arranged in column and rows as in a chessboard. 
The cells have different thickness, however all are chosen to 




Fig. 1. The optical scheme of the nulling interferometer in coax- 
ial configuration that we consider. 




Fig. 2. An example of cellular mirror we introduce in one arm 
of the interferometer. This is just to give an idea : the number of 
cells and their arrangement are arbitrary. 



produce an optical path difference (o.p.d. hereafter) that is an 
odd multiple of Aq/2, where Ao is the central wavelength. Fig. 
[2] gives an example of what a 4x4 cellular mirror could look 
like. Symmetrically, a cellular mirror where cells produce o.p.d. 
which are even multiples of Aq/2 is put into the output beam of 
the second telescope. By design, each cell introduces then, on 
the fraction of the wave at wavelength A® it reflects, a phase shift 
of (2k + l)7r or of 2kn. Hereafter, we call respectively odd and 
even the two cellular mirrors. Of course the number of cells is the 
same on the two mirrors and their geometry is unique (square). 

Note that a Fizeau (or multi-axial) recombination, i.e. done 
at the image plane level and not at the pupil level, could also be 
considered, with the strong a dvantage to keep the design fully 
symm etric. It has been shown (iBuisset et al.ll2 006: Wal lner et al.1 
120041 e.g.) that with this design, nulling is as efficient as in a 
more classical co-axial recombination (Michelson type), pro- 
vided that the star image (i.e. the PSF) is formed on a sin- 
gle mode fiber optics that will make the overlap integral of the 
product between the incident electric fields and the fundamen- 
tal mode distribution of the fiber. In that case, the nulling comes 
from the fact that the amplitude is anti- symmetric with respect 
to the y axis. However, in the following, we will stick to the co- 
axial configuration and especially will show images of the result- 
ing PSF for that type of recombination. The effect of nulling on 
the putative planet detection is seen with a much better contrast, 
most of the stellar light intensity being cancelled. 

The amplitude integrated on the whole pupil area, at the 
output of one chessboard mirror is : 

= Z*Li exp(j0jk) 

where (pk is the phase difference between one cell and a refer- 
ence cell : cpk = 2ndk/A, and dk is the o.p.d. produced by cell k. 
Let's now assume that each o.p.d. is a multiple of Aq/2, where Aq 
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is the central wavelength. The phase then reads : (pk = nrikAo/A. 
If A is within a rather narrow range with respect to Aq, then one 
can write : (pk = nnkil - AA/Aq). 

If now, on one cellular mirror all the o.p.d. are even 
multiple of Ao/2, and on the other mirror, they are odd multiple 
of Aq/2, then the complex amplitudes take respectively the form : 

a+ = Ik=\ expO'7rm£A/i//io)), where rrik are even integers and 

a- = I k=i _ exp(j7rnkAA/ Aq)) 9 where Uk are odd integers. 

When, after reflection on the cellular mirrors, the two beams 
are recombined, then the integrated amplitude a of the wave at 
the output of the recombiner (or, this is equivalent, at the focus 
of the final lens) is : 

a - a+ + a- = ^ =0 exp(j7rmkAA/Ao)) - 
Ik=o exp(j7rn k AA/Ao)) 

Let's note (p = nAA/Ao and let's develop each exponential 
in series of (nk(p) n , so that the two terms of the above expression 
become : 

a+ = I k 1 + j(p xl]k m k-(p 2 x I^k m t ~ J<P 3 x Ik m l + - 
and 

a- = ~ Ik 1 - 70 X Ik^k + 2 xZ^^ + 70 3 x I k n\ + ... 

Clearly, since the number N of cells is the same for the two 
mirrors, the first term of each expression cancels out mutually. 
In order to cancel out the other terms, at least up to a given or- 
der 0, we need to find one set {m^} of even integers and one set 
{tik} of odd integers, such that m° k = ^ k n° k for all o < O. 
Note that several of the m k (respectively of the rik) can be identi- 
cal. This problem is know n as the Prouhet-Tarry-Escott problem 
dBorwein & IngallsTl 994), with here the additional peculiar con- 
dition on the odd/even repartition of the integers. 

The good news is that solutions do exist for any value of O 
! Even more, there are several solutions for a given value of O 
and the number of solutions increases very rapidly with O. For 
instance, let's consider a nulling up to O = 2 : we are looking for 
at least one set of even integers m k and one set of odd integers 
Uk, such that : 

Ik m k = Ik n k and m l = Ik n \ 

One solution (that can be found for instance by trying all possi- 
ble combinations of first odd and even integers) is : 
{m,} = {2,2,2,4},K} = {l,3,3,3} 

One checks that 2 + 2 + 2 + 4=1 + 3 + 3 + 3 = 10 and that 
2 2 + 2 2 + 2 2 + 4 2 = l 2 + 3 2 + 3 2 + 3 2 = 28 

At this point, one remarks that an infinity of solutions can 
be built from a given one : it suffices to add any integer to all 
elements of the two sets. For instance the sets {rrik} = {4, 4, 4, 6} 
9 i n k} = {3, 5, 5, 5}, obtained by adding 2 to each element of the 
solution given as first example, constitute a valid solution. This 
is so because 

Ikirrik + L)° = I k m° k + oLm° k ~ l + fytfrrif 2 + ... + L° , where 
any elementary even sum Yuk (j)L l m^~ l is equal to the equivalent 
odd sum : (°i)L l nl~ l 

From the optical point of view, it's clear that all those solu- 
tions are totally equivalent : this is nothing but just adding every- 
where a same phase shift (or piston). In the following, we will 
use this property of invariance by translation. Hereafter we'll 
consider only solutions where at least one cell has a null thick- 
ness (0) : this means simply that it corresponds to the reference 
level of the mirror surface. 

In the following section, we examine how proper solutions 
can be built in a coherent and systematic way for any value of O. 




Fig. 3. Example of a pattern of cells that is valid for an on-axis 
nulling but that will clearly have drawbacks. The X-Y distribu- 
tion a) will produces aberrations (e.g. tilt) while the b) one is 
more homogeneous and will be preferred. 

2.2. The various solutions 

As we said, there is a rapidly increasing number of solutions 
when O increases. For instance, for O = 3, a systematic search, 
gives 16 solutions when the number of cells per mirror is set 
to 32. However, all those solutions are not equivalent from the 
instrumental point of view. We are looking for solutions which 
are best adapted to our optical problem. One natural condition is 
that the dispersion of the thickness on the surface of the mirror 
is not too large : obviously if the cellular mirror comprises a few 
cells with thickness much larger than the other, then there will 
be both a problem of manufacturing of the mirror (steep edges) 
and of mutual shadowing by the cells. For instance, the solution 
{m k } = {4, 6, 6, 0, 0, 0, 0, 0} , {n k } = {1, 1, 1, 1, 1, 1, 1, 9} 
is valid for (9 = 2, but the odd mirror exhibits a tower of high 
thickness emerging from a flat pond and this is clearly not desir- 
able. 

To avoid large dispersion in thickness, one can put as a con- 
dition that the ordered set of thickness, increases by one unit 
from one cell to the other. In other terms we want to use all 
thickness in a set of incrementing integers. This is a guarantee 
of a rather low variance of the thickness on the surface of the 
mirror. 

2.3. The optimized solution 

We found an elegant and singular solution which satisfies this 
maximum flatness condition and which is based on the Pascal's 
triangle (the rule to derive recursively the binomial coefficients). 
More precisely, for a nulling up to we consider mirrors with 
2° cells, the number of cells of thickness k being the binomial 
coefficient {fy. A clear way to present the solution is simply to 
write the first integers fro to O on one line and the binomial 
coefficients of order O on the next line. 



o.p.d. introduced by cell 





1 


2 


3 




O-l 


Number of cells 


1 


O 


(?) 


(?) 








The first line gives the o.p.d. introduced by the cell (in unit of 
the elementary step Aq/2) and the second line gives the number 
of cells of this kind. For instance, if O = 6 : 



Thickness of cell 





1 


2 


3 


4 


5 


6 


Number of cells 


1 


6 


15 


20 


15 


6 


1 



This means that the even mirror is built as follows : a) it is 
made of 0/2 + 1 (if O is even) or (O + l)/2 (if O is odd) subsets 
of cells; b) all cells in the subset k have an o.p.d. equal to k\ c) 
the number of cells in the subset k is given by the coefficient of 
the binomial, (fy. All even values of k, such that < k < O are 
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thus found at the surface of the even mirror. For instance, 1 cell 
of null o.p.d. , (° 2 ) = 0(0 - l)/2 cells of o.p.d. 2, (° 4 ) cells of 

0. p.d. 4, 1 cell of o.p.d. O if O is even, or O cells of o.p.d. 
O - 1 if O is odd. Symmetrically, the odd mirror is built with the 
same rule, but with k being odd; it exhibits thus : O cell of o.p.d. 

1, (£) cells of o.p.d. 3, (°) cells of o.p.d. 5, 1 cell of o.p.d. O 
if O is odd or O cells of o.p.d. O - 1 if O is even. Note that there 
are 0/2 (if O is even) or (O + l)/2 (if O is odd) subsets of odd 
cells. 

One can demonstrate that this solution do satisfy the condi- 
tion of nulling '.Y^k m °k ~ ^k n l f° r a ^ ° - O ( see Appendix). 

3. X-Y distribution of cells in the pupil : which 
criterion ? 

If the optimum solution provides the distribution in z of the 
cells, it does not say anything on the distribution in x and y, 

1. e. on the surface of the mirrors. Once again, some common 
sense can be exercised to select one with respect to another. For 
instance, if all the large thickness are on one side of a mirror, or 
if there is a continuous gradient of thickness across the surface, 
then this is equivalent to introducing some tilt on the wavefront, 
a bias that is not desirable (see Fig (fright for an example 
of configurations with a same z-distribution but different x-y 
distributions, one being better than the other). A first condition 
that can be stated is that the global tilt be zero. This can be 
expressed by requiring that the first moments of the phase of the 
reflected wavefront is zero : 

2, j i x mi j = and £ Z j j x m Uj = 0, 

where i, j are the coordinates of a cell in a chessboard pattern, 
the origin being at the center of the mirror. This does not seem 
to be a constraint hard to fulfill. To go further, one can ask that 
higher moments of the phase be equal to zero, so that higher 
order aberrations are cancelled : 
Zij i° X mij = and Zij f X m Uj = 

In a final step, all moments, including the cross terms could also 
be considered : 

Z iJ i°xfxm i j = 

Another criterion is to consider the image plane, i.e. the 
plane conjugated to the pupil plane through the Fourier trans- 
form (in an on-axis recombination as above-mentioned). We can 
look for patterns of cells, such that the nulling for the starlight 
extends as far as possible from the centre, for instance along 
the two main axis, with of course the less of attenuation for the 
planet. 

Let's assume that the cells on the even and odd mirrors are 
distributed according to repectively a matrix P and a matrix Q, 
where each element e z -j is the integer determining the o.p.d. of 
the cell (/, j). 

At a wavelength X ± A , the Fourier transform ( hereafter 
noted *{}) of the combined pupil reads : 
f {a+ + a+} - {(exp jPcp - exp jQcp) * cell], 
where cell means a function whose value is 1 inside a square of 
the size of one cell and outside. Note that with this condensed 
notation, exp jPcp means in reality the 2D Dirac's comb with a 
period of one cell width and an amplitude exp jPk,i<P> Using the 
properties of the convolution *, this expression becomes : 

<{a+ + a + ] = f {exp jP<f> - exp jQ(p] x PSFC, 
where PSFC is the point spread function given by a unique 
cell. Since the term PSFC corresponds just to the envelope 
that modulates the resulting PSF, we can restrain the study to 
'{exp jPcp - exp jQ(p}, the term which carries the information on 



the distribution of the cells. Because we are dealing with Dirac's 
combs, the Fourier transforms reduces to a discrete sum : 

r {exp7'P0-exp7'20} = ^k,i (exp j^y-exp jq Ki ) exp ja(kx+ 

ly), 

where pkj and q^i are the elements of matrix P and Q respec- 
tively, and a is a coefficient resulting from the Fourier transform. 
After a Taylor's development, this expression can be written : 

Z*,z IMPkj - qk,i) ~ (p 2 (p 2 kJ ~ q 2 kJ ) ~ j<P 3 (plt ~ q 3 k j) + 
...] exp ja(kx + ly). 

Considering this expression, we see that a sufficient condition to 
get a null contribution at any point of the y axis (x = 0), is to 
have simultaneously : 

Zk(Pk,i - qk,0 = o, I,k(pli - ili) = o> likiplj - vlj) = o, 

etc. for any value of /. 

We remark that this is the same condition, applied here to each 
column, that we found for the whole pupil, on the power of sets 
of odd and even numbers, i.e. Yuk m l - n l f° r a U o < O. 

This means that if we can distribute the cells along the 
columns of the odd and even mirror, so that for each column 
the condition Y,k m l ~ 2fc n °k ^ s sa tisfi e d up to an index O, then 
the amplitude is not only very low at the centre, but also along 
the y axis. Of course, the same result applies along the x axis if 
the condition can be satisfied for all lines. This does not guaran- 
tee that outside those two main axis the nulling will be as deep, 
but by continuity of the PSF, there is a good chance that it re- 
mains fairly good. We will see below that we found indeed such 
distributions of cells that satisfy this condition, both along the 
lines and the columns (see also Appendix D for a different but 
complementary approach). 

3. 1 . The optimized solutions 

We found several ways to geometrically distribute the cells on 
the mirror's surface, following one or the other of the prescrip- 
tions established in the previous section. We describe hereafter 
three of them, each with its advantages and limits. The first one 
corresponds to the first criterion specifying that the first mo- 
ments in the pupil must be zero, while the two other patterns 
satisfy the second criterion on the nulling along the x and y axis 
in the image plane. Other solutions, which are essentially vari- 
ants of those three have been also explored and the section on 
performances will include them. 

3.1.1. First method 

This first method tries to fulfil the condition on minimizing the 
moments on the wavefront after the reflection on the chessboard 
mirror; it is based on a recurrence scheme to build the pair of 
mirrors of order O from the pair of order 0-1. Let's assume that 
three square matrix B, P, Q are defined with dimension 2°~ l x 
2°~ l . The even matrix is P and the odd one is 2, while B is an 
auxiliary matrix. We then build the matrix of dimension 2° x 2° 
as follows : 



Po 


Po-i 


Qo-i - Bo-i 




-Qo-i + Bo-\ 


-Po-i 


Qo 


Qo-i 


P o-i - Bo-i 




-P -i + B -i 


-Qo-i 


Bo 


Bo-i 


-Bo-\ 




-Bo-i 


Bo-i 



As a starting point we choose : P=0, Q=l, B=l. So that at 
the first step : 
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Pi 






Qi 


1 -1 
1 -1 


Bi 


1 -1 
-1 1 



Fig. S] shows an example of the pattern of cells obtained with 
that method and Fig|7]-a displays the resulting PSF with a planet 
10" 6 fainter than the star. 

One can verify that the first moments of the wavefront, up to 
the order O, do satisfy the condition of being null. 

3.1.2. Second method 

The aim of this distribution is to produce an efficient nulling 
along both the x- and y-axis in the image plane. 

Here, we start from the remark that if one writes the first 
2° integers in their binary form, and associate to each integer 
the number of bits set to 1 , the frequency of those numbers are 
distributed according to the binomial coefficients. For instance if 
= 3: 



Integer 





1 


2 


3 


4 


5 


6 


7 


Binary 


000 


001 


010 


011 


100 


101 


110 


111 


#of 1 





1 


1 


2 


1 


2 


2 


3 



The occurence of bits set to 1 is indeed (Tj = ^) = 3 for one or 

two bits set, and Q = = 1 for zero or three bits set. 

If we replicate the elements of this line through all permu- 
tations, we will keep the property of the Pascal's triangle. We 
can then build easily a square matrix with 2° lines of this kind. 
However, we wish to do it so that we also produce columns that 
share the same property. One practical way to obtain this result 
is to write i © j (exclusive OR) at the intersection of line i and 
column j. The number of bits set to 1 on the resulting number 
gives then the o.p.d. of the cell in steps of A /4. An example is 
given below : 





000 


001 


010 


011 


100 


101 


110 


111 


000 





1 


1 


2 


1 


2 


2 


3 


001 


1 





2 


1 


2 


1 


3 


2 


010 


1 


2 





1 


2 


3 


1 


2 


011 


2 


1 


1 





3 


2 


2 


1 


100 


1 


2 


2 


3 





1 


1 


2 


etc . 



















One can check that : a) each line or each column contains 
the right number of cells following the Pascal's triangle rule; b) 
the lines and columns are distributed exactly the same way (the 
matrix is symmetric with respect to the second diagonal). After 



this operation is done, the odd and even numbers are interleaved. 
The odd and even mirrors are then built by simply adding 1 to 
respectively the even cells or the odd cells, as illustrated in the 
following table : 



even mirror 


odd mirror 





2 


2 


2 


2 


2 


2 


4 


1 


1 


1 


3 


1 


3 


3 


3 


2 





2 


2 


2 


2 


4 


2 


1 


1 


3 


1 


3 


1 


3 


3 


2 


2 





2 


2 


4 


2 


2 


1 


3 


1 


1 


3 


3 


1 


3 


2 


2 


2 





4 


2 


2 


2 


3 


1 


1 


1 


3 


3 


3 


1 


2 


2 


2 


4 





2 


2 


2 


1 


3 


3 


3 


1 


1 


1 


3 



We know, thanks to the property of invariance by translation, 
that the condition on the sum of powers is kept. Because of that 
and because of property a) and b), the searched condition on the 
sum of powers along each column or each line is thus satisfied. 

Fig. \5\ shows an example of the pattern of cells obtained with 
that method and Fig(7]-b displays the resulting PSF with a planet 
10" 6 fainter than the star. 



3.1.3. Third method 

The third method aims also at producing an efficient nulling 
along both the x- and y-axis, but is now based on a recurrence 
scheme which is applied to a single matrix Ro (with Rq = 0), as 
follows : 



Ro 



Ro-i Ro-i + (- 
Ro-i-(-l) Ro-i 



For instance : 



Ri 





-1 








1 









Ri 





-1 


1 







1 





2 







-1 


-2 





-1 







-1 


1 






To extract respectively the odd and even matrix, one simply 
add 1 to the even elements or to the odd elements in matrix R, as 
in method 2. One can check that, again, the searched condition 
on the sum of powers along each column or each line is satisfied. 
Fig.[6]shows an example of the pattern of cells obtained with that 
method and Fig|7]-c displays the resulting PSF, assuming a planet 
10" 6 fainter than the star. 

In fact, two versions of this pattern exists : one which is anti- 
symmetric and produces an assymmetric PSF and one which is 
symmetric and produces symmetric PSF. It seems, for practical 
reasons (e.g. alignment of the optical setup) that the second pat- 
tern would be preferable. 
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Fig. 6. Example of pattern of cells in the case of the third 
method, with = 6. The grey level codes the o.p.d. of the cells. 




Fig. 7. Comparaison of the planet image obtained with three dif- 
ferent methods for generating the chessboard mirror pattern. A 
log scale is used. The planet, which has a flux 10" 6 times fainter 
than the star, is the pointlike source at centre; it is surrounded by 
residuals of the starlight. The same cuts have been used, so that 
the difference in intensity seen on the planet is meaningful. 



4. Simulation and predicted performances 

The most important question is now to evaluate if the achro- 
matic phase shifter we propose is indeed efficient, both in terms 
of nulling and of wavelength range. For that purpose We have 
developed a numerical simulation to assess the respective perfor- 
mances of the various configurations we have found. Essentially, 
it consists in : a) building a pair of phase masks for a given type 
and a choice of dimensions ; b) assuming a source exactly on 
axis, to perform a FFT on the sum of the complex amplitudes on 
the two masks, for a given shift of wavelength with respect to 
the central one : this produces the resulting nulled image of the 
star, of course after taking the square of the modulus ; c) to add 
a phase shift, generally n, on one mask, in order to mimic the 
presence of a planet that would be on a constructive interference 
fringe and again to perform the FFT and take the square of the 
modulus. 

In Fig. |7] we compare the different images of the planet - 
assumed to be at a level of 10" 6 of the stellar flux - and residuals 
of the star, obtained fot the typical case (A - A )/A = 0.2. 

The assessment is then done : a) on the quality of the nulling 
versus the wavelength, for instance in computing the residual 
flux in a circular diaphragm centred on the star position; b) on 





Fig. 8. Plot of the theoretical gain in planet to star contrast (log 
scale), as a function of the relative wavelength difference (A - 
A )/A for different x-y configurations of the cells. Top : mirrors 
with 8x8 cells; bottom : mirrors with 64x64 cells. The singular 
peak at the central wavelength corresponds to a perfect nulling 
(gain oo) so that the point translate only the numerical accuracy. 
The different curves correspond to the different patterns of cells. 
The two curves noted 5 and 6 correspond to the two versions of 
the third method (see text). Curve 3 and 4 corresponds to two 
variants of method 2 and curves 1 and 2 to variants of method 1 . 



the magnitude of the unavoidable partial extinction of the planet; 
c) on the wavelength range where a good nulling is obtained. 

The simple metric we have used to compare the different 
configurations is to compute the ratio of the residuals of the 
planet flux to the star flux, each being measured within one res- 
olution element of a single telescope, i.e. an area of A/D around 
the main axis in the image plane. Another possible metric would 
be to estimate a signal to noise ratio on the planet detection, as- 
suming either photon noise (planet + stellar residuals) or a noise 
proportional to the stellar residuals since it seems clear today 
that a lack of stability in the nulling will prob ably b e the main 
source of noise in nulling interferometers iLayl {2004)). 

We present in Fig. [8] the results of the model in terms of 
nulling efficiency vs (A - A )/A for the different configurations 
we have tested and for two numbers of cells in the phase mask. 
The following remarks can be done : 
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First, it is specially rewarding to see that the specifications 
on nulling efficiency that a DARWIN type mission would re- 
quire (a contrast of 10" 6 and in fact a specification of 10" 5 on 
the nulling efficiency), could be easily met with the solution we 
propose : e.g. theoretical contrasts as high as 10 15 are predicted 
on a significant wavelength domain. 

Second, the performance in terms of achromaticity is fair, 
since the level of 10" 6 is reached on the domain [0.6 A Q - 1.25 
A ] with the best configuration : this corresponds to more than 
one complete octave in wavelength. It means that with only two 
devices it is possible to cover the whole wavelength range of 
DARWIN. We note also that there is an asymmetry of the useful 
wavelength domain : the performances are within the specifica- 
tions in a broader domain on the short wavelength side. 

Third, the number of cells is important as shown by compar- 
ing curves at top (8x8) and at bottom (64x64) in FigJSl However, 
increasing this number beyond 64x64 does not bring a very sig- 
nificant gain. 

Fourth, among the different patterns of chessboard, clearly 
the third method presents much better results - at least when 
we consider the raw numbers - since, for instance, it gives 
a contrast 10 6 times better than the next best configuration at 
(A - A )/A = -0.2 . This does not mean that practically it would 
always be the best and a more thorough analysis is needed, tak- 
ing into account deviation from a perfect system, once the device 
is manufactured. 



5. conclusion 

We have presented in this paper a new concept of a quasi- 
achromatic phase shifter, a functionnality which is at the heart 
of any nulling interferometer. We think that it can simplify the 
design and the setup of the interferometer since its relies on 
a simple component of bulk optics symmetrically introduced 
in the two arms of the interferometer. In that sense, it could 
overcome certain limitations of other solutions which use multi- 
components systems and/or complex asymmetric design. The 
heart of the concept is a chessboard mirror with a large number 
of cells (typically 64x64), put in the pupil plane. It is the pe- 
culiar distribution of the different cell thickness which is at the 
origin of the quasi- achromatic capability : it allows to cancel out 
the first terms of the Taylor's development in AA/A of the ampli- 
tude, up to any given order. Among the various distributions that 
satisfy this condition, a very unique one is especially optimum. 
We also show that there is a peculiar x-y distribution of the cells 
which gives very good performances in terms of planet to star 
contrast. Reaching a very deep nulling on a broad wavelength 
range is a mandatory condition for detecting terrestrial planets 
with nulling interferometers and ultimately characterize them, 
so that this solution can be of interest only if it reach the proper 
performance of nulling in a broad enough domain and, indeed, 
we show that the specification of 10" 6 in contrast that a space 
mission as DARWIN is aiming at, can be reached on a complete 
octave in wavelength in the best case of a 64x64 chessboard mir- 
ror. 

Finally, we want to introduce a new terminology. The set of 
equations on which this new concept is based is of diophantine 
type (i.e. polynomial equations where solutions are integers), as 
was a similar set of equations in the question of very deep nulling 
interferometers that one of us is exploring (Rouan 2007) since 
a few years; as a consequence, we propose to call diophantine 
optics this new type of optics where peculiar relations between 
powers of integers are used to obtain very specific effects. 
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Appendix A: The z-distribution based on Pascal's 
triangle 

We want to demonstrate that the condition Yuk m °k ~ ^k n l f° r 
all o < O is f ullfilled with the solution given by the Pascal's 
triangle. We use a recurrence reasoning. The nulling relation we 
want to demonstrate reads : 

zto(- l )^C) = 0forallo - a 

The expression on the left can be rewritten: 
Ito(-l)^- 1 ^) 

We use now the relationship k(fy = O^l) to transform this 
expression in : 

We note that when o ^ the summation can start at k = 1 , so 
that the equation to verify is : 

gc-i)^- 1 ^;*)^ (A.i) 

Now let's assume that the nulling relation holds for p - 1 and 

Let's use the property demonstrated in section 1 that any transla- 
tion of the o.p.d. keeps the relation true, so that, with a translation 

This can be rewritten : 
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This is the same expression as in Equ lA.U where it is just multi- 
plied by O. This demonstrates that if the relation holds for p - 1, 
then it holds for p. It's then sufficient to find a value of p for 
which the equation is true, to prove that it is true for any p. For 

p = 0, the condition becomes : ZfcLoC - !)^) = 0» wnere me left 
part is actually the development of (1 - 1)° and is obviously 0. 



where (p = 7rA. 

The device can be made achromatic up to order O if we si- 
multaneously solve the above system of diophantine equations 

t=i m k = Et=i «* 
ZL «2 = Hi »£ 



Appendix B: Image in the focal plane 

In order to assess the device's performance after insertion of the 
phase grids into the beam, it is convenient to calculate analyti- 
cally the image at the telescope focus. This image is obtained, at 
wavelength A, by squaring the Fourier transform module of the 
pupils. One obtains formally : image(M, v) = module 2 (u, v). We 
have : 

module(w, v) = sinc(^Md) sinc(^vd)|Sum(expUV • expPQ)| . 

(B.l) 

In this expression, u and v are the coordinates on the focal plane, 
d is the size of the mirror, no is the size of the phase grids and 
the definition sinc(x) = sin x/x was adopted. The remaining term 
involves the two square matrices expUV and expPQ. The oper- 
ation • between the two matrices is the elements by elements 
product not the matrices product. The term Sum(A) is just the 
sum of all the elements in a matrix A and, as usual, \z\ = ViE 
z being the complex conjugate. The two above matrices in (IB. lb 
have format no x no and are given by : 



^ m? = Y»nl 



-*=1 k 



jmidd/no ~ j7rvd6/n 



expUV = e 
expPQ = e j* PA e -j™ D + e MQA e +jnuD _ 



(B.2) 
(B.3) 



We have adopted the convention that the exponential of a matrix 
is the matrix of the exponentials. Now, D is the distance between 
the two telescopes and S is a vector of integers varying from 
-(no - 1) to +(fto _ 1) by step of 2. For example if no = 2 we 
have S = [-3,-l,+l,+3]. The symbol is the tensor product, 
accordingly expUV is a square matrix constructed as a tensor 
product of two vectors (sometimes called the outer product or 
the cartesian product). Finally the chromatic term A is: 



1 



^0 



AA = A- A 



o • 



(B.4) 



The above expression is valid for the star of the system. 
For the planet, one must correct imaged, v) with the relative 
monochromatic luminosity at wavelength A and, if one neglect 
the planet tilt on the telescope pupils, add the phase shift on one 
grid. Typically, we set the Q grid to Q + 1, in order to obtain 
maximum planet detectivity through a n phase shift. 

Appendix C: First order: the Prouhet-Pascal 
solution 

At first order approximation, we set D = and S = 0, the cor- 
responding interferometer is an on-axis one, having all cells of 
each phase shift grids on top of the others. In this case, we get 
expUV = 1 and expPQ = e jnPA + e jnQA , leading to the chromatic 
term : Sum = Sum(e jnPA + e jnQA ). As discussed earlier its Taylor 
developpement is given by : 

N N N 

Sum = Yjjnl ~ 4) + M ~ n\) + (j<P) 2 J](m 2 k - n\) + . . . 



k=i 



k=\ 



k=\ 



This system is called a multigrad system and is equivalent to the 

abbreviated notation : {mk,k = l,N} = {n^k - l,Af} or simply 

{mk} = {nk} if there is no ambiguity. Two sets of integers, solu- 
tions of a multigrad system, are called multigrade sets, for ex- 

3 

ample we have {-11, -3, 3,11} = {-9, -7, 7, 9}. The problem of 
solving a multigrad system corresponds to the classical Prouhet- 
Tarry-Escott problem. Our problem is a particular multigrad sys- 
tem since the two sets {mk} and {^} must have the same size and, 
more important, the mk must be even and the nk must be odd. 
The solution of our problem hinges around the propertie that 

if {ak} = {Pk} then {ak,fik + c} °= l {fii,ai + c} for any c. (These 
properties can be easily proved by induction using the fact that 
the highest power coefficient of a polynom in x is invariant under 

the translation x + c.) Starting from {-1, 1} = {0, 0}, we conserve 

2 

parity by setting c = 1 and we get {-1, 1, 1, 1} = {0, 0, 0, 2}. For 
the next step, we regain symmetry by setting c = -I and get now 

{-1, -1, -1, -1, 1, 1, 1, 1} = {-2, 0, 0, 0, 0, 0, 0, 2}. What we use 

in fact is {-1, 1} = {0, 0} and the property that if {m^} = {n^ then 

0+2 

{nk-l,mk,mk,nk+l} = {m^- l,^,^,m^ + 1}, the parity being 
preserved. 



Appendix D: Second order: distributing the cells on 
the pupil 

In the preceding section, we built a nulling interferometer with 
the even cells and the odd cells one on top of the others. The 
practical realization of such an interferometer is conceivable, but 
is far too complicated and will not further be considered here. 
The cells must therefore be distributed spatially on the even and 
odd pupils. 

Take for example of the two 4x4 grids (no = 4). 
The Prouhet-Pascal solution tells us to distribute the two 
sets {-1, -1, -1, -1, 1, 1, 1, 1} and {-2, 0, 0, 0, 0, 0, 0, 2} on the 
pupils, but does not tell us where. To gain some insights on this 
difficult problem, it is necessary to scrutinize the expUV matrix 
a step further and look at its expansion. We have: 

expUV= Yl/p)[UV]*, (D.l) 
^— ' k\ y no' 



k=0 



with (always for the example no = 4 case): 



[UVf = 



(-3u + 3v) k (-u + 3v) k (u + 3v) k 

(-3m + v) k (-u + v)* (u + v) k 

(-3u - v) k (-u - v) k (u - vf 

(-3u-3v) k (-u-3v) k (u-3v) k 



(3m + 3v) k 
(3m + v) k 
(3m - v)* 
(3m - 3v) k 



(C.l) 



(D.2) 

We observe that the even order matrices [UV] 2/? are axisymetric 
while the [UV] 2/?+1 odd ones are anti-axisymetric. Therefore it 
easy to nullify all the contributions from the odd ones in (IB. II) , 
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it suffice to distribute the cells in an axisymetric manner. This 
cancels out the first order term in (ID. 11) and all the odd ones but, 
of course, the even ones are still contributing. Next, we observe 
that the [UV]* matrices are 2-D polynomials in u, v coordinates, 
therefore they will be cancelled out by a 2-D finite differential 
operator. Following these lines we start with the 2x2 matrices 

corresponding to the doubled {-1,1} = {0, 0} solution: 









, Q = 


-l +i 








+i -l 



Indeed, expPQ « P - Q is a differential (gradient) operator. We 

obtain the next matrix by using the {nk - 1, m^, m^, n& + 1} °= 2 
{ntk - l,rik, rik, + 1} property. Since there were still some de- 
grees of freedom left, we found the best arrangement keeping 
the P - Q differential operator property in order to minimize the 
star/planet ratio of the residual light. We got : 



+2 








-2 




+l 


+1 


-l 


-l 














, Q = 


+1 


-l 


+1 


-l 














-l 


+1 


-l 


+1 


-2 








+2 




-l 


-l 


+1 


+1 



The algorithm of the third method proceeds following these lines 
from sl no x no matrices to the next 2no x 2no one. This algorithm 
is very efficient because it produces at last very low star/planet 
ratio. Moreover it possess the desirable property to produces 
axisymmetrical images which helps finding the location of the 
planet. 
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